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' We study the normal modes and the stability of two component condensates in a phase separated 

, regime. In such a regime the system can undergo a quantum phase transition upon the change of 

04 ■ interaction strength between bosons of the same species or the variation of the trap frequencies. In 

^ ' this transition, the distribution of the two components changes from a symmetric to an asymmetric 

5^ , shape. We discuss the nature of the phase transition, the role of the interface tension and the 

■ phase diagram. The symmetric to asymmetric transition is the simplest quantum phase transition 
that one can imagine. We found new branches of normal modes which are specific for trapped 
multicomponent condensates and are analogous to the waves at the interface between two layers of 
immiscible fluids under gravity. At the point of the phase transition the frequencies of those modes 
go to zero and become imaginary which causes an instability. The interface tension shifts the normal 

■ mode frequencies and changes the stability region of the system. 
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I. INTRODUCTION 



■ There is much recent interest in Bose-Einstein condensation (BEC) in trapped gases as well as quantum phase 
transitions (QPT). Examples of the QPT include the Wigner electron solid melting transition, the Mott-Hubbard 
' ' metal-insulator transition, and different magnetic transitions. In physical phenomena involving BECs quantum me- 
^ I chanics play a crucial role. In this paper we propose a new kind of quantum phase transition in phase-separated 
Q . mixtures of BECs. In this transition, the distribution of the two components changes from a symmetric to an asym- 
' metric shape. To explore the transition, we first investigate the stability of the symmetric phase by studying its 
normal modes. We find interface modes that become soft. These unstable modes can be caused by changing the 
I ' trapping potential and are analogous to gravitational waves at the interface between two immiscible fluids. When the 
^ I lowest frequency becomes zero, the instability sets in which determines the stability limit of the symmetric phase. 
' We determine the actual phase boundary by comparing the energy between the symmetric phase and the asymmetric 
. phase and find that the actual phase boundary and the instability boundary is not the same. This suggests that 
' the transition is first order. The system may be a good laboratory to study issues of quantum metastability and 
I . tunnelling. The symmetric to asymmetric transition is the simplest quantum phase transition that one can imagine. 
P\J ' Careful study of this problem should provide us new insight into this burgeoning field of discovery. 

' Systems of multicomponent condensates were first realized by the JILA group in a magnetic trap in ^^Rb and 
subsequently in an optical trap in ^^Na (spinor condensate) [^JJ]. Experimentally at low fields, the spin exchange 
process can occur in an optically trapped condensate, leading to spin domains Q| with metastable behavior [^,^. 
S Binary condensates in two hyperfine levels of *^Rb have been studied 0,0], most notably realizing a system of 

interpenetrating Bose fluids ||], measurements of phase dispersion and a vortex state in a dilute-gas BEC 
(~| Theoretical treatment of such systems began in the context of superfluid helium mixtures ||lO|| a nd spin-polarized 
O ' hydrogen and now has been extended to BEC in trapped gases (l|Jl|,|l|,|l|,|l|,|7|jl8ljl|,^,y| 



^ ' The equilibrium density distributions of segregated mixtures in the absence of gravity have been studied numerically 
for different system parameters. Two types of configurations have been discussed: a symmetric [^,^0[^ configuration, 
for which one component is inside the other one, and an asymmetric one in which the two components occupy the 
left and the right hand side of a sphere [ pO|Jl4|. Several theoretical papers have described collective excitations 
of multicomponent condensates. Busch et al. ]22|] studied collective excitations in the limit of weak interactions 
based on a trial function approach. Collective excitations of binary mixtures in the Thomas- Fermi (TF) limit have 
been investigated by Graham and Walls |^^. Esry and Greene based on the Hartree-Fock and random-phase 
approximation have numerically calculated the low-lying excitations of double BECs in a time-averaged orbiting 
potential trap in which gravity separates the centers of the two components. Gordon and Savage [ p5[ have studied 
excitation spectrum as a function of number of particles. Mazets has performed an analysis of wave dispersion 
on a boundary between two weakly segregated untrapped BECs. 

Under a change of the interaction strength the two component BECs can undergo different types of macroscopic 
quantum phase transitions. For example, an increase in interspecies interactions a\2 results in a transition from a 
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binary mixture to a phase separated state. Such phase separated condensates are analogous to a system of immiscible 
fluids with positive surface tension. Pu and Bigelow p^ ] numerically studied the frequency of collective excitations as 
a function of interspecies interactions ai2 and found appearance of imaginary normal mode frequencies at the phase 
separation transition point. Such imaginary modes describe the nature of system's instability. The trapping potential 
is not crucial for this type of phase transition and similar imaginary modes also exist in the homogeneous (untrapped) 
BECs 

However, in the phase separated regime the trapping potential changes the system's symmetry and results in a 
new type of quantum phase transition. Such transition occurs in trapped condensates upon the change of relative 
intraspecies interaction strength 011/022. When this ratio differs significantly from one the less repulsive component 
is in the middle of the trap and the more repulsive component comprises outer shell. One can make the initial 
configuration unstable by manipulating the ratio 011/022 close to one, e.g., by means of Feschbach resonances or by 
changing the ratio of the trapping frequencies. In this paper we study the collective excitations of the trapped phase 
separated condensates as a function of 011/022 and show that normal modes with imaginary frequencies appear when 
this ratio approaches one. Effect of the interface tension results in a shift of the transition point of the symmetric- 
asymmetric transition in the two component BEG. In general, the asymmetric phase possesses a lower interface 
energy. On the other hand, since the degree of self-repulsion may differ between the two species, the less self-repulsive 
component will prefer to remain where the density is higher, while the other component moves to the low density 
regions outside. This favors the symmetric phase. Depending on the system parameters, one of these two energetic 
considerations will win out. These system parameters can be adjusted by changing the trapping frequencies, the 
relative particle numbers of the two species, and the interaction between the particles with Feschbach resonances. We 
first address the stabihty of the symmetric phase. 



II. BASIC FORMALISM 



Let us consider the two-component BEG in a spherically symmetric trap. The dynamics of the system is described 
by time dependent Gross-Pitaevskii equations 

= -|-A*i -I- Vtr*i + ^^(oiil^iP + ai2|*2n^'i, (1) 
Ot Zm m 



= -:^A*2 + Vt,-^2 + ^^(0221*2!' + ai2|^'in«'2, (2) 
ot Zm m 

where ^'1,2 are the condensate wave functions, Vtr — mujQr'^ /2 is the trapping potential, loq is the trapping frequency, 
r is the radial spherical coordinate, o^ > are s— wave scattering lengths. Here we suppose that trapping frequencies 
and particle masses are the same for both species. We discuss the general case in the last section. 
We shall assume in this paper that the condition 

ai2 - aiia22 > (3) 

is satisfied and, therefore, the condensates are phase-segregated. The first demonstration of a condensate binary 
mixture by Myatt et. al. 0] produced overlapping condensates of the \F — l,mf = —1 > and \F = 2,mf — 2 > spin 
states of ®^Rb. These states, however, possess different magnetic moments and, hence, the condensates experience 
different potentials in a magnetic trap which results in unequal displacement from the trap center by gravity. Later 
JILA experiments were performed on mixtures of 1 > and |2,1 > states These two states have essentially 
identical magnetic moments, and feel identical confining potentials. In ^''Rb, the scattering lengths for 1 > and 
|2,1 > states are known to be in the proportion on : 012 : 022 1.03 : 1 : 0.97, with the average of the three 
being 55(3) A [Q. Hence, the condition (^) is valid for the JILA experiments and the condensates are in a weakly 
segregated phase One should mention that for earlier experiments with the |1, — 1 > and |2, 2 > rubidium states 
the scattering lengths are on : 012 : 022 :: 1.007 : 1 : 1.01 For such states the condition (^) is not fulfilled and 
the mixture exhibits behavior of a miscible system. 

We suppose the condensates can be well described by the TF approximation. In this regime, the phase-segregated 
condensates overlap over the length scale A = H \J a,\2 j ^ a\\ai2 — 1, where f is the healing length [|l^. For the JILA 
experiments on phase-segregated states A « 47^. If, however, the penetration depth A ^ i?, where R is the size of 
the system, the condensates can be approximately treated as nonoverlapping, which we assume to be the case. The 
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effect of the overlapping results in finite surface tension and can be included via boundary conditions at the interface. 
However, in the TF limit such effect is small, we shall discuss it in Sec. IV. If the condensates do not overlap one can 
neglect the last terms in Eqs. (|l|), (||). As a result, the dynamical equations for and \l/2 decouple. However, the 
two condensate components are coupled by boundary conditions at the interface which require continuity of pressure 
and the normal velocity. 

For definiteness we assume that stationary configuration is spherically symmetric with the central core dominated 
by the first component and an outer shell from the second species (see Fig. 1). The stationary density distribution 
"■i = I'I'iP of two components is given by 

= 0<r<R., (4) 



"'"^V ^*<'^<^2, (5) 

where Ga = ATrh^aa/m, Ri — \/2/Xi/ma;Q. The normalization condition J riidV — Ni, where Ni are the numbers of 
condensate particles, determines the chemical potentials fii. The position of the phase boundary i?» is given by the 
condition that pressures exerted by both condensates are equal |29]: 



where we have introduced dimensionless parameters 



/an . M2 

, A=— . (7) 

V 0.22 Ml 




Fig. 1. Symmetric configuration of BEC components. 



The symmetric configuration is favorable when k differs from unity, with the less repulsive component being in the 
middle (an < a22, that is k < 1). At the interface between the two components the stationary condensate densities 
are 

K(Ai2-Mi) (M2-Mi) /o^ 

Gii[l-K) G22(1-k) 

n2/ni = K < 1. 

One can rewrite the time dependent Gross-Pitaevskii equations in a hydrodynamic form which shows an analogy 
between our problem and the motion of two immiscible fluids. In the strong phase-segregated regime the dynamics 
of each components is described by the following hydrodynamic equations 
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_+V.(nV) = 0, (9) 

lmV^ + Vtr~^ ]=AV^+Gn + m^ = fi, (10) 

2 2m y/n ot 

where V is the condensate velocity and $ is the velocity potential, V = V$. Eq. (^ is continuity equation for 
compressible flow, while Eq. ( p^ ) is the analog of the Bernoulli equation with {e + P)/n ~ —^^A^/n + Gn, where 
e is the internal energy density and P is the pressure. The kinetic energy (quantum) pressure is omitted in the TF 
limit and, therefore, e = Gn^/2. As a result, from Eq. ( p^ ) we obtain for the pressure 

2P ~ Gn^ — /in — nVtr — mn-^ — —mnV^. (11) 

The trapping potential plays the role of gravitational potential in hydrodynamics. The linearized hydrodynamic 
equations determine the condensate normal modes. In terms of perturbation in the condensate density n' and the 
velocity potential $ the normal-mode amplitudes satisfy the coupled equations Q 

iun' = V • (Vn') + V • (nV$), (12) 



G 

iuj^ = V • V$ H n' • 



nV — 



(13) 



where n is the static condensate density and V is the static condensate velocity. In our problem V = 0, hence, 
n' « iwTO$ / G and the equation for the perturbation in velocity potential is 

_ 2r^ + (i?2 - r2)A<i> = 0. (14) 
uJq or 

One can seek solution of this equation in the form 

$ = $(r)rz„,(0,0), 

where Yim{9, (p) is the spherical harmonics. Then, using AYim{0, 0) — + l)/?'^, we obtain the following equation 
for the radial function $(r) 



9,9 9, 6?^^ , n 9^ d$ 

dr'' dr 



l{l + l)R'~r' ( ^ + /(/ + l) 



$ = 0. (15) 



Let us derive boundary conditions for $ at the interface. The change in density at the interface is the sum of the 
term from the movement of the interface and that from the perturbation of the density. Let — <;{9, cj), t) be a small 
deviation of the r— coordinate of the interphase from its stationary value R.^,. Then the condensate density at the 
interface is n(i?* + c^, i) « n{R^:, t)+ (,dn{R^,, t)/dr, and, hence, the pressure is given by 

P^^.^^^ + GniR.,tf-^,. (16) 
In terms of the density perturbation n', n{R,:,t) — n{R^) + n'{R^,,t), we obtain 

P « ^a^^^ + Gn{R,)n'{R,,t) - n{R,)mu:lR,^. (17) 
Continuity of the pressure at the interface results in the following boundary condition 

Giini(i?,)n'i(i?»,t) - ni{R^)jnLjlR^(,i = G22n2{R*)n'2{R*,t) - n2{R*)mujlR^<,2- 

Taking into account Gn' = imuj^ and d<;/dt = Vr = d^/dr we finally obtain the following boundary condition for 
the velocity potential at the interface 

LjlR,^-U;Hi^HL0lR,^-KL0^<i>2. (18) 

or or 
Continuity of Vr gives another boundary condition 

fi = ^. (19) 

or or 

Eq. ( p^ ) and the boundary conditions (p^, (19) compose a complete set of equations necessary to determine normal 
modes of the system. 
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III. NORMAL MODES 



For one component condensates the solutions of Eq. (p^ ) which are finite at r = and r — R are 



<I>(r)ocr'Fr-°) 



i?2 



(20) 



where Pn^^'^^^x) are Jacobi polynomials, n is the radial quantum number. The corresponding eigenfrequencies are 
given by 



(21) 



For the two component condensate the normal mode frequencies are different from the general case and new modes 
appear. However, we shall show that frequencies corresponding to n = remain the same as for the one component 
system. Rewriting the function <f> = r^ri{r'^ /R^) in Eq. ( |l5|) results in a hypergeometric equation for ri{^) 



e(e-l)<f + 



1 



(22) 



Hence, general solutions for the inner and outer condensates which are regular at r = and r — R2 respectively have 
the form 



$1 = Cir^F [a, 13,1 + 3/2, , $2 = Car'F U, /?, 1, 1 - 



i?2 

where F is hypergeometric function, Ci, C2 are constants and 
1.3 /,„ . 9 2^2 



Rl 



^+2+\r 



2w2 



Using Eq. (p^), the boundary conditions (|8|), (|T^) and the mathematical identity 



F(a, 13, 7, z) = ^Fia + 1, + 1, 7 + 1, z), 
az 7 



we obtain the equation for the normal mode frequencies 



— = (!-«) 



[/(/ + 3/2) + (/ - cj2/cj2)Axsi(u;,x)] [{I - w2/cj2)a;s2(w,x) - I] 
[l(l + 3/2){K-l)+x{l-ujyuj^) {kXsi{uj,x) + (/ + 3/2)s2(t^, x))] ' 



where 



, , i^(a + l,/3+l,; + 5/2,Ax) , , F (a + 1,/3 + 1, 2, 1 - x) 
^H'^' ^) = FT — ^ , , 0/0 X \ ' *2(w, x) = 



i^(a,/?,; + 3/2,Aa;) 



F(a,/?,l,l-x) 



i?2- 



(23) 



(24) 



One of the solutions of Eq. (24) is 



Iu!q, which coincide with those for one component condensate. For this 



solution the components oscillate in-phase and $1 — k^2 oc r^Yi,n{0 , 4>) . Another exact solution is lu"^ = 5wq with 
$1 cx 1 — 5r2/3i?2j $2 oc A(l — 5r2/3i?2) which corresponds to Z = 0, n 1. For this solution the components oscillate 
out-of-phase if -\/3/5i?i < i?* < ■\/3/5i?2 and in-phase otherwise. 



A. Numerical results 

We solve Eq. ( |2^ ) numerically and find normal mode frequencies ui of the two component condensate as a function 
of the parameter k = ^ 011/022 for different fixed ratios R^jR^- The ratio i?*/i?2 can be directly measured exper- 
imentally. Fig. 2 shows the normal mode frequencies for I — At fixed I there is an infinite number of branches 
which correspond to different radial quantum numbers n. We plot only the lowest modes with n = 0, 1, 2. In Fig. 3 
the normal modes are estimated for 1 = 1. In the limiting case an = 022 (k 1) the condensates behave as an one 
component system and the condensate modes coincide with Stringari's result (p^). Fig. 4 shows the low frequency 
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modes that become imaginary at k > 1. These modes are peculiar for two component systems and are analogous 
to the waves at the interface between two layers of immiscible fluids under gravity | |3^ . In trapped condensates the 
gradient of trap potential plays the role of gravitational field. As soon as k becomes greater than 1, which means 
77.2 > ni at the interface, the system becomes unstable the same way as two immiscible fluids in gravitational field 
when the more dense layer is on the top. One should mention that normal modes with imaginary frequency also 
appear at the transition point from a binary mixture into a phase separated state which occurs at ai2 = ^aiia22 
| p9| , p3|j27| . In that case, the fastest decaying mode occurs at a finite wavevector, resulting in a quantum spinodal 
transition | |3^ . 







1=0 


R*/R2=0.75^^ 




0.25 


n=1 





0.0 0.5 1.0 1.5 2.0 2.5 

K 



Fig. 2. Normal mode frequencies for / — and n = 1, 2 as a function oi k — aii/a22. The mode with n — 2 depends on 
position of the interface, the mode frequency is estimated for _R*/i?2 = 0.25, 0.5, 0.75. For k = 1 the mode frequencies coincide 
with those for the one component condensate. 
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1=1 








n=1 


p /p _Q'75 

R^/R^-O.S 

R*/R2-0.25 


n=0 





0.0 0.5 1.0 1.5 2.0 2.5 

KT 

Fig. 3. Normal mode frequencies for I = 1 and n = 0, 1,2. 




Fig. 4. Low frequency modes as a function of k = •\/aii/a22 for different I = 1,2,3. Tlie position of the interface is 
-R, = K^jl- Frequencies become imaginary at «: > 1. 



B. Low frequency modes 

Here we study the low frequency modes in detail. In the region |1 — k| <C 1 the mode frequencies are small: \uj\ <C wo- 
In this limit one can take w = in the right side of Eq. and put k, A w 1 in the multiple. As a result, we obtain 
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cu^^u;l{l-K)f{l,x), (25) 

where 

_ I {I + 3/2 + xsi{0,x)) {xs2{0,x) - 1) 
^^'""^ x[si(0,:e) + (? + 3/2)s2(0,x)] ■ 

Eq. ( p5[ ) describes behavior of the low frequency modes in the region close to the point of instability k = 1. In this 
region uj cx — n and becomes imaginary when k > 1. In Fig. 5 we plot / as a function of a; = i?*/i?2 for different I. 
Imaginary part of frequencies is greater for larger I and decreases with increasing R^/ R2, that is bigger inner droplets 
are more stable. Figs. 6, 7 show the radial distribution of the density perturbation for the low frequency modes 
with different I. In our estimates we choose the radius of the inner droplet i?* = i?2/2. The density perturbations 
are normalized so that 713(^2) — 1- For the low frequency modes the two condensates oscillate out-of-phase. Fig. 
6 corresponds to k = 0.9 (stable region), while for Fig. 7 k = 1.1 (unstable regime). One can see that the normal 
mode profile undergoes no changes at the point of phase transition, although the mode frequencies become imaginary 
at this point. However, due to the relation n' m iojm^/G, the phase shift between oscillations of n' and <i> changes 
at the instability point from tt/2 to 0. Modes with small I are delocalized and the whole condensate is involved in 
oscillations. With increasing /, however, the modes become localized near the interface and in the limit / — > 00 they 
are similar to gravitational waves at the surface of deep water with the dispersion relation uj oc Vk. 




) ~i — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' 

0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 5. Dependence of the factor / in Eq. ( psj ) on the position of the interface i?«/i?2 for modes with different I = 1,2, 3. 
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Fig. 6. Perturbation in the condensate density of the two components n'i(r) and n'2{r) for the low frequency modes with 
different I = 1, 2, 3. The position of the interface is R* = R2/2 (dash hne) and n = 0.9 (stable region). 



8 
6 
A 
2 
O 
-2 



1=3 \ 


R^/R2=0.5 




K=l.l 


1=2 V 








1=1-^ 














1=1 








1=2 




^ 1=3 



0.0 0.2 O.A 0.6 0.8 1.0 
r/R2 

Fig. 7. The same as in Fig. 6, but for k = 1.1 (unstable region). 

IV. EFFECT OF INTERFACE TENSION 

Interface tension is small in the TF limit and results in corrections of the order of ^/i?2 to the normal mode 
frequencies. However, for the low frequency modes in the vicinity k 1 the effect of the interface tension is substantial 
because the mode frequencies themselves are close to zero. The surface tension a modifies the boundary condition 
for the pressure so that the pressure difference at the interface is equal to the surface tension pressure Pi — -P2 = 
(7(l/ri + l/?'2), where ri, r2 are the principle radii of curvature. The stationary position of the interface is determined 
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by the equation 



2 2 

which results in the following ratio of densities at the interface 



n2 
ni 



4cr 



The interface shape oscillates during the condensate motion which produces oscillations of the surface pressure. For 
interface shape close to spherical [p3 



1 1 

n r2 



2 



2£ _ Ae. 

m 



(26) 



where Ag^^ is the angular part of Laplace's operator. Also one should take into account that cr itself depends on local 
density and varies during the interface motion. We assume a oc 
n{R^ +<;,t)^ n{R^,) + n'{R^,t)+ (^dn{R^)/dr, we obtain 



Giini 



^i^nin2 Then, taking into account 



a = ao [1 



Giini 



1 + 2kp 



2(1 + Kcff) 



m 



G22«2 



(icj$2 - OJqR*<;) 



Kcff 



2(1 + 



where ctq is the tension for the stationary configuration. As a result, the boundary condition ( [Tq ) contains extra terms 
and becomes 



2mCTo 
-R^Giifii 



1 + 2Keff 



mn2 



2(1 + Koff) 

2m(To 



9 9 

LJoR^— 

or 



^(2 



9$i 

9r 



i?*G22f^2 



Kcff 



2(1 + Kcff) 



9„ 9$2 9, 

wgi?,— u;^$2 

or 



(27) 



while Eq. (p^) remains the same. For the radial part of the velocity potential we obtain 



1 - 



CTO 



1 + 2Kcff 



= Koff ( 1 



2(1 + Kcff) 



ujiR, 



R,P 



dr 

2 + Kcff 
2(l + Koff) 



cro(;-i)(; + 2)a$i 



9r 



or 



(28) 



where P is the pressure at the interface. The modified boundary condition results in the following equation for the 
low frequency modes 



1 - 



3(70 , ao(/-l)(/ + 2) 



2R^P 



- Kcff / 



(29) 



where ni is the density at the interface. The interface tension shifts the frequencies of the lowest modes and narrows the 
stability region. For / = 1 the inner droplet moves as a whole without changing its shape. However, the displacement 
of the droplet into the less dense region decreases the interface energy which is proportional to n^/^. This is the origin 
for the contribution — 3(To/i?*P to the mode frequency. The mode frequencies now become imaginary at 



K > 1 — 



2(To(^- l)(/ + 2) 
mniWQPf 



(30) 



As K, increases the mode with I = 1 becomes imaginary first which determines the system's stability limit. When 
the surface tension is absent, the higher I is the more unstable the mode is. In the present case, the two component 
condensate is locally unstable when 



ail ^ ^ ctq 

022 



P,P 



(31) 
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V. GLOBAL STABILITY: ENERGY OF ASYMMETRIC PHASE 



So far we have discussed the local stabihty of the symmetric phase. We next turn our attention to its global 
stabihty by comparing the energies of the symmetric phase and the asymmetric phase. Let us consider an asymmetric 
state in which component 1 sits on the top of the component 2. The transition occurs when k is close to 1 which we 
assume. Hence, the interface is approximately flat (see Fig. 8). 



symmetric 



asymmetric 





Fig. 8. Symmetric and Asymmetric phases of two component BECs. 

In the Thomas- Fermi approximation, the density distribution is still given by Eqs. (Q) and (|^), but the interface 
position is now different. The TF density distribution suggests that the outer edge of the condensates has a spherical 
shape and, therefore, in the asymmetric configuration the condensates have equal chemical potentials /ii^ — /i2a 
Carrying out the integration, we find that the energy of the asymmetric state is given by 



Ea = ^{Nl 



N-, 



1 



1 

G22 



4 
35 



d 
6R 



10i?5 21i?7 



35G: 



Esa, 



22 



(32) 



where z — d is the position of the interface, Ega is the interface energy and R = y^2^/mujQ is the radius of the sphere. 
The radius R (and hence fi) and the position d arc determined by the number of particles 



^1 



N2 = 



From Eqs. (p|), (M), we obtain 



TTflR^ 


4- 


15d 


5^3 3d^ ' 


15Gii 






2¥ 


TTflR^ 


44 


15d 


5d^ 3d^ ' 


15G22 






" "R3" ^ 2^ 


iiiVi + 


G22 


N2 - 


AnmujQR^ 
15 ' 



which determines R. The position of the interface can be found from 

^ Ibd ^ 5d^ 3d^ _ SGiiTVi 



(33) 



(34) 



(35) 



(36) 



i?3 2R^ G22N2 + GiiiVi ■ 

The energy of a spherically symmetric configuration when the first droplet is located inside and the second outside is 
given by 

Rl 



-En = 



^lA^l ^2N2 



Y Ml 


Ai2 ^ 


Rl 


^ mcj2 / I 1 \ 


AGii 


G22 J 




2 \G22 Gil / 
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35G22m5/2cj5 



E 



sO) 



(37) 



where Esq is the interface energy in the symmetric state. 

When the system is close to the point of phase transition the interface energies of two configurations are given by 
(for details see Sec. VI) 



%^hRli,^l^ I , ^ / Rl 



(38) 



Esa — 



5/2 



(39) 



In general, these expressions need to be numerically evaluated. They are analytically tractable in the limit when one 
of the components is particularly small, to which we turn our attention next. 



A. Energy for iVi < N2 



If the number of atoms of the first species is much less than N2 then R^ R, d/ R ^ 1 — (4A^i/5A^2)"'^^^ and from 
Eqs. (H), (|9|) we find 



E 



sO 



-y/ 012/^^^11022 - 1, 



(40) 



imSxI'-^Eso-^Eso. 



The interface energy of the asymmetric state is negligible. Close to the point of phase transition the energy 
difference between the asymmetric and symmetric configuration is given by 



E„ 



n ^2 4V6e / , . 

(1 - k) — Yai2/Voiia22 



R 



(41) 



where ^ = h/^/2mjI. 

In the absence of the interface tension (f 0) the asymmetric configuration has greater energy no matter what is 
the value of k in comparison with 1. However, the interface tension makes the asymmetric phase preferable at some 
AC and results in symmetric-asymmetric phase transition. The condition Ea = Eq describes a line of global instability 
of the system 



On 
022 



1/2 



(012/^011022 - l)^^'*- 



(42) 



The soft mode condition gives the line of local stability 



On 

022 



1 



3-R* 



(012/^011022 - 1)"^/^. 



(43) 



In Fig. 9 we plot the phase diagram that shows different stability regions of two component condensates. In 
estimates we take S./R* = 0.01. When ai2 < ^011022 the homogeneous binary mixture is a stable state. Otherwise 
the two components are phase separated. In the later case for small ratio 011/022 the state (1,2) with the first 
component being inside and the second outside is the only stable configuration. If the ratio 011/022 increases the 
configuration becomes globally unstable when we cross the left (solid) curve. However, the system is locally stable since 
the configuration corresponds to a local minima of energy. Further increase of 011/022 crosses the line of local stability 
(another solid line) and the system undergoes a phase transition into a new stable asymmetric state. If initially the 
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system is prepared in the (2,1) state then with decreasing the ratio 011/022 the configuration first becomes globally 
unstable when we cross the right dotted line and locally unstable when 011/022 is close enough to 1. 
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Fig. 9. Phase diagram of the two component condensate in coordinates displaying relative interaction strength between 
bosons. 

VI. ENERGY AS A FUNCTION OF DROPLET DISPLACEMENT 

The question about the system's global stability is related to the energy of the state. To further explore the 
metastability of the system here we estimate the energy of the system for a displaced inner droplet. We assume the 
inner droplet has a spherical shape with a radius tq and the center of the sphere is displaced a distance zq from the 
trap center. The radius tq depends on zq. The condensate densities are given by the same TF expressions (^), (|^), 
however, the position of the interface is now different. We assume the number of particles of each species are fixed: 



An 
SGii 



10 



(44) 



47r 1 
3 G22 



3 -jmujlrl mu^r^ 
M2ro + ^^ + ^ 



(45) 



For fixed numbers of particles, the chemical potentials are functions of the positions of the droplets and are not 
constants. In the TF limit we find for the energy of the system 



E 



Gil 



JJ2_ 
G22 



5 3 

35G22m5/2tj5 



1 

G22 



1 

Gil 



(46) 



where Es is the interface energy. The energy (]^^ depends on four parameters: /ii, ^21 zq, rg. However, only 
two parameters are independent because there are two restrictions imposed by particle conservation (^), ( |45| ). For 
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simplicity let us consider the path for which fX2 = const. Then only one parameter is independent. If we choose rp as 
an independent variable then the energy as a function of rp can be expressed in an explicit form 

^ , 27r(^io - ti2)^Rt ^ 
E = const H — o h 

(g22 Gil ) 

where /X2 = ma;gi?|/2, /xio is the chemical potential at zq = and i?, is the droplet radius at zq — 0. The dependence 
of the droplet radius on displacement zq is described by Eq. . 



mi 



5R. 



2' 



50i?|r3 



25i?| 



175i?| 



(47) 



A. Interface energy 

Here we estimate the interface energy as a function of the droplet displacement zq. According to Ref. |l^ the 
interface tension is given by 



(48) 



where represents the single condensate coherence length = h/\/2miGiini, the pressure P w Giinf/2 and the 
condensate densities nt are estimated near the interface. The interface tension a depends on the local density and, 
hence, changes along the interface. To calculate a we assume that the system is close to the point of phase transition, 
that is Gil ~ G22 — G, Ri ~ R2 — R, ^J.l ~ iJ-2 — ^J■ and near the interface ni « 77.2 = n. As a result. 



2?iVGn3/2 



and we find for the surface energy 
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(49) 



(50) 



B. Energy for small inner droplet 



In the limit tq 
displacement zq: i 
given by 



R {Ni <C A^2) Eq. ( |45| ) results in a simple relation between the droplet radius tq and its 
^ A'2^2/(a*2 — 'Ti'^o-^o/2)- Then near the point of phase transition the energy of the system is 



where 



E « const + Nifi 



(1 



z2 



5/6' 



(51) 
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\J ai2/%/aiia22 
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2Pi?, 



The first term in Eq. ( |5l| ) corresponds to energy contribution due to difference between an and 022, it increases with 
Zq when K < 1. The second term is the contribution from the interface energy, it decreases with increasing zq because 
surface tension drops when the droplet moves in the region with smaller density. In Fig. 10 we plot the energy (js^) as 
a function of displacement of the inner droplet zq. The energy is shown for h — 0.3 and different k. Far from the point 
of phase transition the energy monotonically increases with displacement of the droplet, in this case the first term 
dominates. When k goes closer to 1 the first term in Eq. (j5^) becomes small and the interface contribution changes 
the energy behavior. The energy becomes nonmonotonic function of displacement. The value of k when the final and 
initial energies are equal is the onset of global instability. However, the system is still locally stable since position 
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of the droplet at the center corresponds to a local minima of energy. In this region an energy barrier prevents the 
droplet to move outside. Further increase of k results in disappearance of the energy barrier and initial configuration 
becomes both locally and globally unstable. 



0.10 




Fig. 10. Energy of the system as a function of displacement of the inner droplet for fixed b — 0.3 and different k. 

Near the point of local instability the energy barrier is AE sa 18Nifi{Ks — K,)^/5b, where Ks = 1 — 56/6. The barrier 
becomes of the order of fi when — k sa ^/5b/l8Nl. If iVi ^ 10^, b ^ 0.1 we obtain Ks — i i - 10"'', while 0.1. 

VII. NORMAL MODES FOR DIFFERENT TRAP FREQUENCIES AND PARTICLE MASSES OF THE 

TWO SPECIES 



In this section we generalize our results to the case uji UJ2, 'nii 7^ ^712 which is useful for future experiments where 
the trap frequencies and particles can differ for each component. For simplicity we omit the interface tension and 
consider spherical trapping potentials. The generalized boundary condition for the velocity potential reads 



miLOiR, 



dr 



miuj $1 — Km2UJ2R*- 



9£2 

dr 



Km2LU $2, 



(52) 



where k = ^ G11/G22 = \J a\\m,2l a22'm\. The boundary condition for d^i/dr remains the same as Eq. (p^. Now 
the phase transition is controlled by a new dimensionless parameter 



K. = 



2 ir^ 'ill 2 I 

m2L02 (jii UJ2 ail 



rriiuj'f V G22 ra^'^iJl V «22 ' 
Instead of (|2^), the solutions for the inner and outer condensates have the form 

r.2 



$1 ^ Cir'F ai, 131,1 + 2,12, 



Rl 



$2 = G2r'i^ a2,/32,l,l 



(53) 



where i?f = ijii/niiuof [i = 1,2) and 
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1 + 



2^2 



,3 L , 9 2cj2 
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Using Eq. ( Js^ ) and the boundary conditions we obtain an equation for the normal mode frequencies 

2 ujI{1 - a) [1(1 + 3/2) + {I -uj'^/ujI)Xxsi{uj,x)] [{I -uj'^/uj'^)xS2iu},x)~ I] 

" mi - ijyLul)Xxsiiu;,x)ijf/Lul + {l + 3/2) {l{kujl/Lol - I) + [l - lj"^ /ujI)xs2{lo,x))] ' ^^^^ 

where 

^ miUjliJL2 Rl , . + l,/3i + lJ + 5/2,Ax) , , F (aa + 1, /52 + 1, 2, 1 - a;) 

rn2t^2Mi ^2 F (q;i,/3i, / + 3/2, Ax) F (q;2, /32, 1, 1 - a;) 

In terms of n, and A the position of the interface is given by 



-R* — i?2 



/ 1-kA 
A(1-k) 



In the region |1 — k| ^ 1 one can take oj = in the right side of Eq. (54) and put k,A « 1 in the multiple. As a 
result, we find 



,2 



ujl{l~k)f{l,x), (55) 



U^ + 3/2 + xsi(0,x))(a;s2(0,a;)- 1) 
■'^'^^ [xsi(0,a:)u;f/w| + G + 3/2)(c^2/^2_^_^^^2(0,x))] ^ 



where 



Eq. ( p5D is a generalization of Eq. (|25|) for the low frequency modes. The point of instability is now determined by the 
condition k = 1, or aiiTO2'^| = a22m\uj\. k can be changed not only by changing the relative interaction strengths, 
but also by changing the ratio of the trapping frequencies. Furthermore it is proportional to a higher (4th) power of 
the ratio a;i/w2 as compared to 011/022 and thus will be very sensitive to the frequency change. Near the instability 
point uj oc \/l — k and becomes imaginary when k > 1. 

In summary, we study the normal modes and the stability of the phase separated two component BECs and predict 
a new symmetric-asymmetric transition in such a system. We show that the normal mode frequencies differ from 
those for the one component condensate and new branches appear. We found normal modes with low frequencies 
which are analogous to gravitational waves at the interface between two immiscible fluids. Under the change of 
the relative interaction strength 011/022 or the ratio 101/102 the frequencies of those modes go to zero and become 
imaginary. Imaginary frequencies mean instability and the system undergoes a quantum phase transition into a new 
stable configuration. Such transition can be observed experimentally by manipulating the interaction strength via 
Feschbach resonances or by changing the trap frequencies. Interface tension shifts the normal mode frequencies and 
results in stability of the asymmetric state in some parameter region. Another interesting possibility is to study 
tunneling by turning on a microwave field which allows direct exchange of the two components and has been done in 
many JILA experiments with two hyperfine states of ®''Rb. If we suddenly switch the species, one can study a decay 
of new unstable configuration or tunneling directly, instead of trying to tune the ratio 011/022. 

We are grateful to A. Fetter for valuable remarks. This work was supported by NASA, Grant No. NAG8-1427. 
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